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Periodically kicked turbulence is theoretically analyzed 
within a mean field theory. For large enough kicking strength 
A and kicking frequency / the Reynolds number grows expo- 
nentially and then runs into some saturation. The saturation 
level Re sa can be calculated analytically; different regimes 
can be observed. For large enough Re we find Re" at oc Af, 
but intermittency can modify this scaling law. We suggest an 
experimental realization of periodically kicked turbulence to 
study the different regimes we theoretically predict and thus 
to better understand the effect of forcing on fully developed 
turbulence. 



Periodically driven flow is ubiquitous. Faraday's ex- 
periment Q is an idealized version thereof, more relevant 
examples are the earth's atmosphere, driven by periodical 
heating of the sun, or the blood flow in veins, driven by 
the beating heart. Another example is the gas flow inside 
a sonoluminescing bubble which is periodically kicked by 
the collapsing bubble wall || . Another example of peri- 
odically kicked turbulence is the numerical realization of 
homogeneous shear flow J3| where periodical remeshing 
is necessary. 

In this paper we set up a mean field theory for peri- 
odically kicked flow, based on the mean field theory for 
decaying turbulence Q which was able to describe the 
experimentally measured energy decay in turbulent liq- 
uid helium flow with fixed external length scale || . The 
goal of this paper is to theoretically understand the dif- 
ferent flow regimes which are to be expected, to explore 
the effect of intermittency corrections on these regimes, 
and to ultimately initiate experiments. 

Another motivation for the paper is to study the effect 
of a specified type of forcing on turbulence. In most the- 
oretical studies on turbulence a Gaussian random noise, 
acting on the largest length scales, is assumed. Only re- 
cently experimentalists started to systematically vary the 
type of forcing || f?J . This work is a further step towards 
the analysis of a more specific type of forcing. 

To define the model, we have to (i) calculate the energy 
input during the kick and (ii) know how the energy is 
dissipated in the time At between successive kicks. 

(i) Kick: As an illustration, consider plane shear flow 
in the 1-direction; the flow is sheared in the 3-direction. 
The width of the channel is L, the velocity of the upper 
plate is U, the lower one is at rest. The average energy 
dissipation rate can be calculated to be (see e.g. H) 
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where (). denotes the average over the x-y-plane and v 
is the viscosity. For laminar flow the first term in the 
bracket is zero and the second one is —vU/L. In turbu- 
lent flow in the middle of the channel the second term on 
the rhs will hardly contribute for large enough Reynolds 
numbers. The first term represents the total turbulent 
flow energy E, order of magnitude wise. Therefore, in 
general, 
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E {t) = - e (t) ~ -E{t) + jpv. 
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Imagine now a short intense kick of time Atkick <C At 
on the flow by rapidly moving the upper plate with U. 
After this kick the initial energy E increases according 
to (H) 
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where we have assumed Atkick <C L/U. Assume isotropic 
turbulence in the flow center and define a Reynolds num- 
ber Q 
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Then eq. (||) translates to 
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di- 



\AtkickUjL < 1 and 



with the 
mensionless kicking strength A 
the "laminar" Reynolds number Rei am — \AtkickV 2 jv. 
We choose this name because for very small Ren we have 
Rei = Rei am . For very large Reg 3> Rei am we have 
Re\ — (1 + A)ReQ. Shear flow and the "derivation" of 
eq. (^]) are only thought of as a motivation; there will 
be many other experimental situations where the energy 
input roughly corresponds to a law of type (^|) . 

(ii) Decay: In ref. (J] we calculated how the turbulent 
activity decays within a time t for flow with fixed external 
length scale L. The calculation was based on Effinger and 



1 Note that this is not the standard definition and gives lower 
values for the laminar-turbulent transition than what one is 
used to. 
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Grossmann's variable range mean field theory of turbu- 
lence |j] in which viscous subrange and inertial subrange 
can be treated equally well. The result of ref. j| is that 
for given initial Reynolds number Rei the time depen- 
dence of Re(t) (defined as in eq. (f|)) is determined by 
the inverse function of 



t(Re) 



[F(Re) - F(Re l )} 
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where r = L 2 jv and F(Re) is given by 

F(Re) = ^ {-7 + Vl 2 + Re 2 } 
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c e ,oo = (6/6) 3 / 2 (in this theory) is the dimensionless 
energy dissipation in the large Reynolds number limit, 
7 = 9/ce.oo, and b is the Kolmogorov constant M] which 
is the only free parameter in the theory of ref. From 
experiment jll[] b — 6.0, a value which we take in all cal- 
culations here. Consequently, c ej00 = 1.0 and 7 = 9.0. 
Rather than the Reynolds number eq. ( 4 ) one could also 
give the Taylor- Reynolds number Re\ \ i \ 
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Eqs. (||) - (0) with t — At (the time between successive 
kicks) fully define the present model. The two main phys- 
ical parameters in the model are the kicking strength A 
and the kicking frequency / = 1/At. The third physical 
parameter is Rei am , the minimal Reynolds number after 
a kick. We pick Rei am = 1 throughout. 

Figure [I] shows Re{t) for fixed A and three different 
kicking frequencies / for the initial Reynolds number 
Rev = Reiam = 1- During each cycle there is a kick 
Reo — > Re\ according to eq. (Js[) and a subsequent decay 
according to eq. (^). Overall, there is growth up to some 
saturation level Re sat (A, /), achieved after t sat (A, /). In 
this saturation state energy input and loss through de- 
cay in At balance and the degree of excitation fluctuates 
between a lower saturation level i?e ; sat and an upper sat- 
uration level 



Re s u at = fleryi + 2A+ {Re lam / Ref at ) 2 . (9) 

The (lower) level of saturation Ref at (A,f) is given 
through the implicit equation 
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FIG. 1. Time series Re(t) for fixed A = 0.1 and three 



kicking frequencies 
bottom to top. The 



are i?ef. 
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f = 10/r, J = 100/r, 'and J = 500/r, 
corresponding (upper) saturation levels 
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a^g-log plot of Re\ a * as a func- 
d eking strength A. Two regimes 
Re\ a ^ < 7 a laminar regime in 
the kick always brings back the 



level of excitatQn |oKt 
For large -Re*"* > y^m nave 'a tut 
with i?e* at cc /. The transition from the laminar to the 
turbulence regimefjakes place around -| 

= jL t/x 
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For large Re^ at , Ref at ^> 7 one has the explicit result 



Similarily, these two regimes are also seen in fig. || where 
we plotted Ref a * vs A. In the turbulent regime for large 
Re\ a * it is Re{^ oc A. 

We now come to the important question of how in- 
termittency effects JlJ] change the exponents calculated 
within this mean field theory. In refs. fl3|| intermittency 
effects have been included into the mean field theory 
of ref. H on a phenomenological basis. One possibility 
for their effect is that the dimensionless energy dissipa- 
tion rate c e becomes slightly Reynolds number dependent 
oven in the large Reynolds number limit jl3|], c e oc Re~ K , 
vith k = (9/8)<5C2/(l + 3<5C2/8). Here, 6£ 2 « 0.03 is the 
ixperimentally found deviation from the classical scaling 
- exponent £2 — 2/3 of the second order velocity structure 
function. The consequences of this small (n w 0.03) scal- 
ing correction can straightforwardly be embodied in the 
present mean field approach to periodically kicked tur- 
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FIG. 3. 



Re s u at and Ref at for / = 10, 10 2 , 10 3 
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FIG. 2. Saturation level Ref at (lower curve of pair) 
and -ReJ a * (upper curve of pair) as a function of_^>for,. 
three different kicking strengths A = 0.01 (solid), <A' =J$.\ 
(dashed), and3tQ 0^5 (dashed-dotted), bottorn-t«'tqp^*l3elow, 
Re a u at ~ 7 = 9.0 the [excited state is laminA^abcam Rel at jxS§ 
it is turbulent. 
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to top. The increasing difference between i?ej and'i?e; at 
the right edge of the figure has its origin in the breakdown i 
requirement A <Z 1. 



rather than ell-.^y)/ Equation @ may offer a new and 
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We expect the totQl energy to builcFup again and again 
over time scales larger than At a^^ho^f'^ suddenly 
drop because a energy pulse is traveling cfownscale. Such 
behavior has been observed in numerical simulations of 
periodically remeshed homogeneous shear flow || and in 
simulations of periodically kicked shell modells of turbu- 
lence [results to be publihsed]. As based on a mean field 
theory the model of this paper is only applicable to the 
mean energy and not to these fluctuations. 

We now come back to an experimental realization. It 
will be easier to perform experiments in a closed system 
rather than in a channel flow. A particular suited exper- 
imental setup for periodically kicked turbulence would 
be the flow in a cylinder between two counterrotating 
disks ||, Q. Also, Rayleigh-Benard convection may be 
well suited. Here, as an example we take Taylor-Couette 
flow If the radii of the inner and outer cylinder 

are similar, the energy input will still roughly follow 
eq. (||). Take water as fluid which has v = 10~ 6 m 2 /s 
and take L = 1cm. Then r = 100s. Realistically 
achievable kick-strengths would be A = At^ickU / {2L) ~ 
0.1s ■ 0.1m/s/ (2 ■ 1cm) = 0.5, e.g., the whole range A < 1 



in which_the theory is applicable^^Por this value the'sc'al- 
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stinguish between eqs. ( pd| ) aitd'( 1^ 
one r^s to ^o up 4e-freqtTencies of around p&z which 
should bE"achievable.__At_ these, rekjtiwly" large kicking 
frequencies measur ements can only reveal instantaneoi 
values. To get st atements on the averaged quantitk 
dealt with in this paper, ensemble averaging is necessary. 
This is best donje by repeatedly probing the flowiat some 
fixe d phase after llg; respective ki£k. Averaging_c|ver the 
results at phase + will give i?e* at ,. averaging over the 
results at phase At~ will give iief ,twio 

We suggest to perform a periodically kicked turbulence 
experiment and to measure Re sat as a function of both A 
and /. To our knowledge it would be one of the first ones 
with some active control on the type of forcing. Immedi- 
ate questions to ask are: Does the level of saturation for 
large Re indeed only depend on the product Af as sug- 
gested by eqs. ([ll]) and (j^) or do boundary effects carry 
on into the strongly turbulent central regime and cause 
a more subtle relation? If so, an application of the so 
commonly used volume forcing for turbulence becomes 
more questionable. Do the (scaling) relations between 
the quantities introduced in this paper, e.g., Ref°^(f,A) 
offer a new way to measure intermittency effects? What 
modifications arise if forcing and decay do not decouple 
as assumed in this simple model? 
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